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(57) Abstract 

!ilmuliu«) CMP results to cwrespondinTcm^^ic^^^^ modeling system opiim Z6S it, modeling pM,n,etc» by comparing *e 

pa«mete« .uch tha. the^STctt ^3t a^^So^a id« ^"^.^ei^J^ r'ji'° capable of opOn^W^ <5S» p4e« 
system leverages Wstoricrafcj data to gene^^e^m^olaSon fomule^^^ ' '^^"'^WS^r. The modeling 

formulas enables the modeu3y„em to simlte CMli"Su^ ItX^ Z IpS'ZTr'*"- ^ 
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METHOD AND APPARATUS FOR MODELINr a 
CHEMICAL MECHANICAL POLltmNGPR^^^^^^ 

TBCHNfrAr FrFTnoFTRFiMVF>m^^f 
The present invention relates generally to chemical mechanical polishing (CMP) system, 
S thatp„ce.woHcpiecessuchassemicon.uctorwafe„. Mo. panicul^ .he eseTL^^^^^^^^^^ 
relate, to a modeling syste. chat simulates CMP .suits and facili.tes :he optimi ^ 
prpcoss p„3 associated with the CMP system in accordance with a peci J C 
between global unifotmity and feature-scale planarity of a workpiece. 

BACKQROTrND OF THF nviyp^rrrj.^ 

' methodfT of semiconductor wafe. has become the prefe.ed 

method tor planara^ng dielectric and other material layers at various stages of integrate, circui 
rab„cat.n. ^^^^g the CMP process._awork^^ 
™ay be covered with one or ^or^^s^^^^^^^^ 

Po.yu.thanepad,Theslur«yusedfor most CMP .y^emsisawater^sedcomli^^^^ 
suspended colloidal silica particles. "ion navmg 

re„,o« acros. «.e waft, .H. U» „ ^^^^^^^ 

a„ a„ e„„.a.a. ,„ ,o,„, 

P^v..,„, ^^^^ P^b 

or PC example, a ngid upper pa. ,e„d. ,o pla„„i.a *e .„fa«Tf I 

P..v,d. eener eloW »,fo„^,, beeaa^e U can epafo™ ,o ,„e ove.„ »h.pe Jc„„„„Tr, ^ 

o;tr.:r.r'' """" -« — t 

op,,™.zauo„ ™, have o„ re..„„.^,e p,a.,„,, „r„„ „„,p,ece. Sue. ,y.Sr2 
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effectively strike an optimized balance betw . PCT/uS98/2423i 

mode/mg systems «r«.- • ^^^'^ P^^^^ss parameters Funh ™P'ncai data 

- capable oZ ""^ ^ ~ CMP syZ S . 

Pbleofprocessmgempiricaland/orsimulatednr. "^^^^^ing systems are 

Accordingly, it is an 
°'P'ovid,„, ^ J .1 ~ «.« .he CMP . 
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to thereby minimize simulation errors. 

Another advantage is that the present invention provides a CMP modeling system that 
employs interpolation techniques to effectively simulate CMP results for which little or no 
empirical data exists. 

An additional advantage of the CMP modeling system is that simulation data collected for 
an existing CMP system may be used to assist in the design of a new CMP system having 
different physical and process characteristics than the existing CMP system. 

BRIEF PFSCRTPTTOM OF THE PR A wr^f q^ 
A more complete understanding of the present invention may be derived by referring to 
the detailed description and claims when considered in connection with the Figures, where like 
reference numbers refer to similar elements throughout the Figures, and: 

FIG. 1 is a schematic block diagram of a chemical mechanical polishing (CMP) modeling 
system in accordance with the present invention in an exemplary operating environment that 
includes a CMP system; 

FIG. 2 is a flow diagram of a CMP simulation process performed by the CMP modeling 
system; 

FIG. 3 is an exemplary wafer-scale simulation result produced by the CMP modeling 
system; 

FIG. 4 is an exemplary feature-scale simulation result produced by the CMP modeling 
20 system; 

FIG. 5 is a flow diagram of a feature-scale simulation process that may be performed by 
the CMP modeling system; 

FIG. 6 is a schematic rendition of a modeled die pattem and a corresponding deformation 
model associated with a polishing element; 
25 FIG. 7 is a flow diagram of a CMP process parameter optimization process performed by 

the CMP modeling system; 

FIG. 8 is a flow diagram of a model validation process performed by the CMP modeling 
system; and 

FIG. 9 is a How diagrajn of an interpolation fomiula generation process performed by the 
30 CMP modeling system. 
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Referrme .o FIG. I, a chemical mecMcd p„fohi„g (CMP) modeling ,„ ^ 

acc.^ a p.fe.e<. em.odimen. o^e p..e„. i„ve„,io„ be us., m ^rl" 
v„^..s.„.eCMPs..em A"houg. mo«i„g .,.«m ,0., CMP ./J 
uauzed .„ .He co„,ex. of of wo^^cce. sys.cm. eomp.Ub,e fo, 

were. .„ de^Hted he„i„ ,he of convenience „d 
.pP"o.««d ,ba. ^ pre^en, invenUon is no. limited „ pa^cula, CMP sy«em or any .peeific 
^ Of wo^Mcce. F.«he™ore. CMP syatem, « generally we,, known in d« aemiconduefor 
fabnca„on mdus,n, and „„, „o, be deacribed in de«i, hc^in «,eep, where necesaary tor an 
understanding of ihe present invention. 

Modeiing system 'Ogenerally includes at leaatacomputer 14. a nserinterface ,6 and an 
output dev,cc „. ,„ ,Ke prefaced embodiment. compu»r >. is con«gmed a. a conv^tionai 
pe-^n^ computer having a processor 20 and a memory 22. in addition to dte ptocesaea 
descrtbcd berem. computer ,4 may perform any number of c„„ve„,i„„ai «„„,io„a unrelated to 
■tlE»«ion and computer i 4 ma, include addidona. h^dware eomponentalSlhS^Jn- 
mFIC Computer I4may be alternately configure! aa pan ofamainftame comparing aya,«n 
par, Of a network environment, or integral to CMP system 12. Indeed, compute, ,4 may be" 
reabzed m a variety of fonns so ,ong as it includes a sufficient amo«n, of computing power and 
memory capacity to support the requirements of the p^aent invention 

Processor 20 is prefentbly confgured to eany out anumber of processes (descdbed below) 
employed by modeling sysrem ,0, Such proeessea may be ,eall»d by software progtamming 
tnatrucons s,o.,d within memory 22 or within a separate memory eiemen. as Jated with 

compu.erM.F„r,he3akeofco„ve„ie„ce,HG.Idepic,saCMPslm„.a.ionprooesa24.amodeI 
vahdauan process 26, a CMP optintization process 28. «mam,!M«H&^^ , 

2S .■.,«,S,,«,..and a featutcscale simulation process 3, aa Z^T^S^^ ' 
processor 20. Processor 20 (,nd/or computer ,4, is p„,brably configured to intent wi«, user 
■merlace I output device 1 8, and CMP system 12 (via, e.g., a CMP contmller 32, 

Memory 22 cooperates wi,h processor 20 in a conventional mam,e, to store, update and 
prov,d. modeling. CMP, and user-de^ned data to »,d fiom processor 20. Although memory 22 
may be arranged in any suitable manner, FIG. . depict memory 22 organized into a number of 
d,s nc, daub..,es for the sake of clarity. r„ partieuto. modeling system lOpteferably inciudes 
leas, . CMP process p_ database 34. a modeling paran«^da!ite^„ 
m«rpol.,,on formula da.a.base3sJ. user-defined CMP ■^^iS-SS^in^:^^^'^!^, 
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T'"" database 42, and a simuJated CMP . '"CT/uS98/2«3, 
r. ^'"^^aWe medium cn^h « more of the 

syslem 12 may com™ ■ moaelmg sys,n„ |o 

" ''^ — --en, of semiconductor iayer rh " --«^u«me„t systems. SystemI 

•0 scale rfnr ^^"'""^ °P«^«« to S'mulate CMP 

" ' -'''"«tes) and a micro T '^^"^^ °" ''o^h a wafer 

Modeling system in • ^ «na a micro-feature scale rfor « , 

recine • " """'"^ of optimizing the CMP . ^ ^ ""''"^"^ ^^^^^'"^e^)- 

recpe. ,n accordance with a user-defin.H • • "^'''^ Pa^meters. . 

Drofilp "ser-aefined initial die oafffm process 

profile, and/or a user-defined vveigj,tin« of ^^^-defined initial JocaJ feature 

"-fonnity. Modeling ..^tem iO jreslT"^^ ^'--^^V and ^ 

~inaccordanee.i.histol^^^^^^^^ 

-"•a.e CMP processes for which no empi W ^ "^'"'"^ ^"""'^^ -d to 

data e.sts or to simulate the performance of 

Refernng now to FIG 2 CMP • 
embodiment, an operator is n 

through user interface J6 r »u 
5 °' ^"'"es generated by 
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modeling system 10, or accept one or more default values presented by modeling system 10. 
Process 24 preferably includes a task 52, during which modeling system 10 obtains a 
. plurality of CMP process paratneters associated with a CMP procedure to be performed upon a 
workpiece. In the context of this description, a "CMP process parameter** is any quantity, 
5 characteristic, dimension, or other variable that may have an effect upon the outcome of the CMP 
process. For example, such CMP process parameters may be related to the desired polishing 
recipe, e.g., polis h time, speed of thep olishin^ elem ent and associated accelerationj^amp time, 
workpiece carrier speed and acceleration ramp time, caiTier down force and application ramp 
time, carrier sweep range, and carrier sweep speed, A different polishing recipe may also be 

10 specified for distinct polishing stages of the CMP process. In addition to polishing recipe 
information, the CMP process parameters may relate to the particular CMP system 12 utilized 
during the CMP process, e.^,, dimensions and/or hardness of the polishing pad, dimensions of 
the workpiece to be processed, characteristics of the cairier sweep range and pivot point, and the 
like. It should be noted that any niunber of CMP process parameters may be obtained during task 

IS 52 and that modeling system 10 may provide one or more default CMP process parameters or 
facilitate the selection of one or more preexisting CMP system configurations during task 52. 
The CMP process parameters may be stored in database 34 (see FIG. 1) for future use. 

CMP simulation process 24 also involves a task 54, during which a plurality of initial 
modeling parameters arc obtained by modeling system 10. Once obtained, the modeling 

20 parameters may be stored in database 36 and subsequently accessed by modeling system 10. 
These modeling parameters are utilized by modeling system 10 to compute a simulated CMP 
result associated with, inter alia, the CMP process parameters. In the preferred embodiment, 
each modeling parameter is expressed in terms of a quadratic equation having a number of 
user-definable equation coefficients. Thus, task 54 may obtain or calculate the initial modeling 

25 parameters by receiving a number of such equation coefficients from the user. In the context of 
the present invention, each of the modeling parameters may be a function of one or more of the 
CMP process parameters, including: the polishing table speed, the workpiece carrier speed, the 
workpiece cairier down force, the polishing element composition or construction, the slurry 
characteristics, and the particular film layers being processed, each of which is associated with 

30 the CMP procedure to be simulated. 

Although the present invention may incorpomte any number of suitable modeling 
techniques, the preferred embodiment utilizes methodologies derived from the Preston equation 
for material removal: 
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Material Removal Rate (R) = kPS, 
where k (Preston's coefficient) is a factor representing chemical effects, P is the pressure 
imparted by the wafer onto the polishing pad, and S is the relative speed between a point on the 
wafer and the polishing pod. In other words, the product PS represents the mechanical effects 
associated with the CMP process. Although this model is relatively simple, the value of it has 
a significant amount of inherent uncertainty because the chemical effects represented by k 
include the chemical reactions between the polishing slurry and the wafer in addition to the 
availability of the polishing slurry at the surface of the wafer. Consequently, k can be affected 
by the composition of the polishing slurry, the composition of the wafer, the characteristics of 
Ihe polishing pad, the wafer carrier speed, and tlie speed of the polishing table. Indeed, although 
k may be held constant to simplify the simulation routine, it may actually vary across the surface 
of the wafer. 

Although P may also varj' across the surface of the wafer, it is typically better defined than 
k. For example, tlie average value for P can be computed by dividing the carrier down force by 
the surface ai ea of the wafer. Additional factors such as the curvature of the wafer backing film 
and positions of vacuum holes (employed by the wafer carrier to secure the wafer) may be 
analyzed to effectively estimate how P is distributed across the surface of the wafer. For the 
exemplary CMP application described herein, P varies quadratically with the radius of the wafer. 

'I'hc computation of the relative speed betwee n e^^f Tp oint^ the wafe r and the underlying 
pad is accomplished through a straightforward applicationofkmematics. Accordingly, the value 
of 5 can be predicted with near certainty. Because S does not affect * or P, its effect on tlie CMP 
process can be predicted with a high degree of confidence. 

In the preferred embodiment, k is one of the modeling parameters; changing k results in 
a corresponding scaling of the material removal rate associated with the simulated CMP process. 
In addition to k, modeling system 10 preferably employs a number of additional modeling 
parameters to specify the pressure distribution across the surface of the wafer. As described 
above, the actual average pressure across the wafer can be easily calculated if the carrier down 
force and surface area of the wafer are known. Thus, modeling system 10 uses tlie additional 
modeling parameters (up to five in the prcfened embodiment) to approximate the actual pressure 
distribudon across the wafer and to simulate the corresponding material removal characteristics 
associated with the CMP procedure. It should be noted that modeling system 10 may be 
configured to hold k equal to zero when the current interrogation or sampling point is located 
over a groove or gap in the simulation polishing element or when the simulation wafer is 
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overhanging the edge of the simulation polishing element This additional consideration may 
be desirable to compensate for those sampling points where the polishing element does not or 
cannot contact the wafer 

Modeling system 10 may be configured to provide one or more default modeling 
parameters for a given CMP simulation; the user may have the option to alter such default 
modeling parameters as desired. Similarly, modeling system 10 may provide a number of default 
coefficients for a given modeling parameter to simplify the amount of custom data entry required 
during task 54. As described in more detail below, task 54 may simply obtain a set of modeling 
parameters that have been optimized for a particular CMP system, a particular wafer size, and/or 
a particular wafer composition. 

Still referring to FIG. 2, CMP simulation process 24 may also perform a task 56, during 
which a number of simulation control parameters are obtained by modeling system 10. The 
simulation control paranieters are related to the manner in which the CMP simulation is carried 
out. For example, modeling system 1 0 is preferably configured such that a simulated CMP result 
includes data associated with a plurality of discrete interrogation points on a given workplace. 
The use of distinct sampling points is desirable to enable precise comparison of the simulated 
CMP results to empirical CMP results measured at the same points. Accordingly, the simulation 
control parameters may be related to the location and number of sampling points on the wafer, 
the coordinate system used for an output plot, and the time periods associated with repeated 
simulations for common sampling points. Of course, any number of additional settings or 
variables related to simulation protocols may be obtained during task 56. 

A task 58 is preferably performed to obtain user-defined CMP data associated with the 
CMP procedure to be analyzed. The data obtained during task 58 is preferably stored in database 
40 (see FIG. 1) for use by modeling system 10. For purposes of the present invention, 
"user-defined CMP data" means theoretical, desired, or actual characteristics of the original or 
processed wafer that may have an effect on the CMP simulation. For example, the suggested 
CMP recipe generated by modeling system 10 may be dependent upon an initial feature-scale 
pattern associated with the workpiece, an initial local film thickness profile associated with the 
workpiece, a desired level of global wafer uniformity, or a desired level of wafer planarization. 
In addhion, the simulation procedure conducted by modeling system 10 may be dependent upon 
an indicator of the relative importance of global wafer uniformity versus local die planarization 
for the workpiece. 

During task 58; data indicative of an initial die pattern and/or an initial film thickness 

8 



Wo 99/25520 

PCT/US98/24231 

profile may be obtained in any suitable format. For example, a given pattern may be fonncd 
from a plurality of arrays each defmed by a plurality of nodes; each array may be defined by the 
nuniber of features contained therein, the length of its head section, the length of its tail section, 
the mmiber of "hills" and the height of the hills; the lengUi of the spaces or "valleys" between the 
hills, and other descriptive elements. In practice, the user may input physical characteristics of 
the pattern features, numerical descriptors, or a graphical rendering of the pattern to convey the 
initial die pattern to modeling system 1 0. Similarly, the initial thickness profile may be rendered 
i n any suitable manner, e.g. , on a point-by-point basis. 

The indicator of uniformity veVsJis planarity may be expressed as^^p. a percentage, a 
.•scaled number, a graphical representation, or in any suitable manner. For instance, an evenly 
weighted indicator may cause modeling system 1 0 to optimize both uniformity and planarily and 
generate a particular CMP recipe to reflect the even balance. Hovvever, an indicator that favors 
uniformity may cause modeling system 1 0 to produce an entirely differem CMP recipe that is " 
intended to increase the uniformity of the wafer at the expense of feature-scale planarization. 

It should be noted that all of the above tasks 52. 54. 56. and 58 need not be performed 
during CMP simulation process 24 and that the above tasks 52. 54. 56. and 58 may be petformed 
in a different order than that described herein. In addition, the relevant data need not be 
organized, arranged, or obtained in the specific manner described above. For example, the CMP 
process parameters, the initial feature-scale pattern, and the indicator of global wafer uniformity 
versus local die planarization may all be categorized as CMP data. Likewise, the modeling 
parameters and simulation control parameters may all be categorized as modeling or simulation 
data, in addition, the data and informadon set forth above may be received by computer 14. 
memoty 22. or processor 20 (or other components of computer 14) in any suitable manner and 
such data and infonnation may eventually be stored in memory 22 or routed to processor 20 for 
further manipulation in accordance with the various processes carried out by modeling system 
10. 

After the appropriate data is received by modeling system 10. a task 60 causes modeling 
.system 1 0 to perform an appropriate modeling routine to obtain a simulated CMP result for tlie 
given workpiecc. The simulated CMP result is eventually stored in database 44 (see FIG. 1); 
this data may be used during subsequent processing by modeling system 10. In the preferred 
embodiment, the simulated CMP result is associated with the CMP process parameters, the initial 
modeling parameters, and at least one element of the user-defined CMP data. In other words, 
the simulated CMP result generated by modeling system 10 may be sensitive to changes in any 

9 
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of the user-defined parameters. 

As described above, one exemplary governing equation relates the film thickness at a 
number of discrete sampling^points^o^an^ associated with the CMP procedure. 

In this example, let (x^.y,) be the coordinates of sampling points on the surface of the wafer 
(where / = I, 2, 3, . , . » jV), and let be the film thiclcness at the particular sampling point. Tlien, 
the current film thickness for eacff sampling point is governed by the following differential 
equation: 

where: 

/=L2,3 N; 

k(t) = Preston's CoefiFicicnt; 

P = Pressure at the sampling point; 

== Velocity of the pad'^under the sampling point; and "5^ ^ 
V„. = Velocity of the wafer at the sampling point. 

In practice, modeling system 10 performs a number of numerical computations to solve 
the above differential equation based on Preston's relationship (and possibly other equations); 
the resulting solution is related to the theoretical material removal rate associated with the 
specified CMPl5rocedure. The theoretical pressure distribution derived from the modeling 
parameters is applied during the simulation to predict how much material will be removed at the 
specific sampling points on the surface of the wafer. It should be appreciated that the rate of 
material removaL which is dependent upon the amount offeree applied by the polishing pad to 
the wafer, may differ from die to die and within each die. Modeling system 10 preferably 
analyzes the amount offeree applied to the wafer on a localized scale and determines the erosion 
of the wafer in response to the localized distribution offeree. 

During task 60, modeling system 10 preferably obtains at least a wafer-scale simulation 
result and a feature-scale simulation result for the current workpiece. In an exemplary 
embodiment of the present invention, the wafer-scale simulation result includes a simulated film 
thickness profile and the feature-scale simulation result includes a simulated local pattern profile. 
The film thickness profile may be utilized to derive global wafer uniformity information, and the 
local pattern profile may be used to derive local die planarization information. The .global wafer 

10 
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uniformity information is related to the simulated film thickness calculated at the various 
sampling points, while the planarity information is related to tlie localized flatness of the pattern 
at particular die locations. Feature-scale simulation process 31, which may be performed by 
modeling system 10 to obtain the feature-scale simulation result, is described in more detail 
5 below. The simulated CMP result may contain additional information related to the wafer 
characteristics and the simulated CMP result may be formatted or expressed in any suitable 
manner. 

After the modeling routine is completed for the wafer (or after several modeling runs are 
completed over a particular time period), a task 62 causes modeling system 1 0 to produce a 

10 suitable output for review by the user. As described above in comieciion with FIG. 1, output 
device 18 may be a display terminal, a printer, a plotter, or any suitable device capable of 
conveying the appropriate simulation results to the user. In the preferred embodiment, task 62 
produces an output indicative of the wafer-scale simulation result and/or the feature-scale 
simulation result. FIG, 3 is an exemplary wafer-scale simulation result (e.g., a film thickness 

15 profile) produced by modeling system 10 and displayed on a conventional computer terminal. 
Each plot 70 represents the film thickness at certain locations on the surface of the wafer, e.g., 
sample points taken along the diameter of the wafer. Different plots 70 for the same simulation 
may be associated with the condition of the wafer at different processing times. For example, 
a plot 72 may represent the film thickness at time while a plot 74 may represent the film 

20 thickness at a time / + 

FIG. 4 is an exemplary feature-scale simulation result (eg., a local pattern profile) 
produced by modeling system 10. This particular simulation result includes an initial user- 
defined feature-scale pattern 76, e.g., the original local die pattern obtained during task 58 (see 
FIG. 2) or a derivative tliereof. The simulation result may also include one or more plots 78 - 

25 representing various stages in the simulated CMP procedure. Such simulation results may be 
utilized to determine the effect that CMP process parameters have on the local planarization of 
the wafer. 

After task 62 produces one or more outputs related to the modeling results, CMP 
simulation process 24 ends. It should be appreciated that process 24 may continue with any 
30 number of additional tasks and that process 24 may be incorporated into one or more 
comprehensive processes utilized by modeling system 10 or CMP system 12. 

FIG. 5 is a fiow chart depicting feature-scale simulation process 31 . Process 31 may be 
performed by modeling system 10 during CMP simulation process 24. Process 31 preferably 

1! 
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begins with a task 150» which causes modeling system 10 to pbtain the initial feature-scale 
pattern associated with the workpiece. As described above in connection with task S8 (see FIG. 
2). the initial feature-scale pattern may be stored in memory 22 for subsequent access during 
process 31. Referring to FIG. 6, an exemplary feature-scale pattern 170 is illustrated. Pattern 

5 170 is shown as if the wafer is face-down in a position ready for processing by a polishing 
element. As shown, feature-scale pattern 170 is preferably represented by a plurality of nodes 
172 connected by line segments for purposes of simulating the CMP procedure in accordance 
with the present invention. The spacing between nodes 172 may be selected in any suitable 
manner, ^.g., to provide a sufRciently accurate model of pattern 170, 

1 0 Referring back to FIG. 5, a task 1 52 may be performed to acquire a deformation model of 

a polishing element, e,g,, the polishing pad used by CMP system 12. The deformation model 
may be a user-defined model, a default model employed by modeling system 10, or a 
combination of user-defined and default characteristics. An exemplary deformation model 174 
is shown in FIG. 6 in a deformed state that simulates the characteristics of an actual polishing 

15 clement used by CMP system 12. In .the preferred embodiment, deformation model 174 is 
defined at nodes 176 that correspond (relative to the^honzojatal axis) to nodes 172. Further, 
deformation model 174 may assume that adjacent nodes 176 are connected by line segments. 

At any given time during the CMP simulation, feature-scale pattern 170 is considered to 
be rigid and the polishing pad is considered to be deformable in accordance with deformation 

20 model 174. Deformation model 174 is preferably defined by a plurality of primary 
force/displacement elements 1 78 that are capable of compressing in response to a load. In the 
exemplary embodiment shown in FIG. 6» each primary force/displacement element 178 is 
associated with one of nodes 176. Each primary force/displacement element 178 may be 
"linked" or otherwise associated with at least one additional primary force/displacement element 

25 178. In the preferred embodiment, primary force/displacement elements 178 are "coupled" lo 
each of their adjacent primary force/displacement elements 178 via secondary force/displacement 
elements 180. Any number of modeling parameters may be employed to characterize primary 
and secondary force/displacement elements 178, 1 80. Displacement elements 178,. 180 may be 
considered to function like springs for purposes of analysis. 

30 After task 1 52 acquires pad deformation model 1 74, a task 1 54 is preferably performed to 

initialize a displacement 1 82 of the wafer relative to the polishing pad. In practice, displacement 
1 82 may relate to an amount of downward travel associated with the wafer carrier during the 
CMP procedure. The initial displacement may be set at any suitable value. In response to tliis 

12 
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displacement, a task 156 causes modeling system 10 to determine (by simulation) the 
deformation of the polishing element. During task 156, modeling system 10 may suitably 
generate a simulated contact profile and/or a localized force profile associated with deformation 
model 1 74 in relation to the current state of feature-scale pattern 170. The contact profile may 
5 identify which nodes 172 of pattern 170 are in contact with corresponding nodes 176 of 
deformation model 174, As shown in FIG. 6, portions of the simulated polishing element may 
not be in contact with the simulated wafer die pattern, and secondary force/deflection elements 
1 80 may limit extension of some primary force/deilection elements 1 78 that would otherwise 
extend to contact the simulated die pattern. The force profile may then identify the amount of 

10 localized force imposed upon nodes 172 by nodes 176 by analyzing the appropriate 
force/deflection elements 178, 180 with respect to displacement 182. 

The^wMiiact profile is preferably computed by utilizing static force equilibrium e^ations 
at nodesA76yi^or any given node, the downward displacement of its primary force/deflection 
element induces a conesponding upward force. The relative displacement of this node with 

1 5 respect to its adjacent nodes induces additional forces (either upward or downward, depending 
on the relative positions). Modeling system 10 preferably assumes that the summation of the 
nodal forces for any given node is equal to zero. When that limitation is imposed mathematically 
for a node, an equation in terms of the displacement of the given node and the displacement of 
the nodes adjacent to the ^ivcn node is formed, ^^o^inin^uch^^ for each of the nodes 

20 1 76 forms a sec of algebr aic^tfations in terms of the displacements of nodes 176. In that set of 
equations, the force-balance equation for any node that has previously been determined to be in 
contEp?with"the ^ is replaced with an equation holding that the position of the node is equal 
to thq^ verticaljrosition of the wafer at that point, 

Afterthe system of^equatio ns is solved, the polishing pad displacements are examined. If 

25 a node position is determined to be higher than that of the corresponding die pattern at that point, 
then the node position is automatically set to be equal to the die pattern position and is identified 
as a "possible contact node," The system of equations is then re-solved, replacing the force 
balance equation for each "possible contact node" with a specified displacement, which is- the 
position of the wafer surface at tliat node. Similarly, if a node, previously thought to be a 

30 ''possible contact node" is, based on the force balance, determined not to be in contact, then it is 
removed from the list of "possible contact nodes" and its corresponding equation is replaced 
with the force balance equation. The solution of the equations and subsequent determination of 
"possible contact nodes" is repeated until no changes occur between iterations. At that point, the 
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"possible contact nodes" arc determined to be in contact with the wafer and are identified as 
"contact nodes.'* The force imposed on the wafer by each contact node is computed using the 
force relationships associated with primary and secondary force/deflection elements 178. 180. 
After task 156 generates an appropriate force profile for the current displacement 1 82, a 
5 query task 158 is performed. Query task 158 preferably comparesjhe sum of the nnrfps for/^^g . 
contained in the localized force profile to the current (simulated) down force associated with the 
pkrtjcultti Lareffc-efTEe wafer being analyzed. This down force may be a user-defined quantity^or 
it may be derived from a known carrier down force associated with the entire wafer and the 
known area of the wafer under analysis. If query task 158 determines that the sum of the node 

10- forces does not substantially equal the local down force (or, alternatively, if the difference 
between the sum and the down force does not fall within a desired tolerance), then a task 1 60 is 
performed to adjust displacement 1 82 by a suitable amount. The adjustment of displacement 1 60 
causes a change in deformation model 1 74, the current contact profile, and the current localized 
force profile. Accordingly, task 156 and query task 158 are preferably repeated for the updated 

IS displacement 182. 

Tasks 1 56, 158, and 160 preferably form a processing loop that causes displacement 1 82 
to be adjusted in a suitable manner until query task 158 determines that the sum of the node 
forces is approximately equal to the local down force applied to the wafer. If these forces 
substantially balance, then a task 1 52 is performed to cause modeling system 10 to simulate the 

20 erosion of the wafer for a given time period* In the preferred embodiment, erosion is simulated 
by adjusting the positions of the appropriate nodes 172 by an amount that may be dependent 
upon the current CMP parameters, the slurry composition, tlie wafer composition, and the like. 
Thus, for a given time soon after CMP processing begins, the local die pattern may closely 
resemble the original die pattern (see FIG. 4). At a later time, modeling system 10 may leverage 

25 historical simulation data to enable simulation of erosion over time. For example, a subsequent 
iteration of feature-scale process 3 1 may utilize a partially planarized version of the die pattern 
rather than the initial die pattern described above in connection with task 150. 

FIG. 7 depicts optimization process 28 that may be performed by modeling system 10 to 
generate a preferred set of CMP process parameters for use during an actual CMP procedure. 

30 The CMP process parameters are preferably optimized to produce a CMP polish recipe in 
response to an intended CMP result. Process 28 may be carried out in conjunction with CMP 
simulation process 24 or as a separate and distinct process. 

Optimization process 28 preferably begins with a task 86, which causes modeling system 
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10 to retrieve a current simulated CMP result from, e.g.^ database 44 in memory 22 (see FIG. 
1). Following task 86, a task 88 is performed to retrieve the optimization parameters for process 
28. These optimization parameters may be retrieved from database 40. In the context of this 
description, an optimization parameter may be any characteristic, quantity, or feature associated 
5 with an "ideal" wafer as processed by the specific CMP system 12. The preferred embodiment 
of modeling system 10 may include any number of optimization parameters related to: the initial 
feature-scale pattern; the initial film thickness profile; the relative importance of wafer 
uniformity versus wafer planarization; or the intended planarization or uniformity result. These 
parameters may be user-defined. 

10 Following task 88, a task 90 is performed to cause modeling system 10 to compare the 

characteristics of the current simulated CMP result to the corresponding characteristics indicated 
by the optimization parameters. For example, in the context of feature-scale planariiy, task 90 
may analyze an error between the intended die pattern and the simulated die pattern contained 
in the current simulated CMP result. Similarly, with respect to wafer-scale uniformity, task 90 

1 5 may analyze the difference between the intended film thickness and the simulated thickness 
measurements. 

A query task 92 is preferably performed to determine whether the error between the 
simulated CMP result and the intended CMP result is substantially minimized. Optimization 
process 28 may utilize any number of teclmiques to analyze the simulation error, e.g., curve 

20 fitting, least-squares, averaging, and the like. If query task 92 determines that the simulation 
error is substantially minimized, then the current simulation is considered acceptable and query 
task 92 may prompt a task 94, which saves and displays the current CMP process parameters. 
Following task 94, optimization process 28 ends and the operator may apply the optimized CMP 
process parameters to the actual CMP procedure. If query task 92 finds that the difference 

25 between the simulated and intended CMP results is not minimized, then a task 96 may be 
performed to cause modeling system 10 to adjust at least one CMP process parameter. For 
example, task 96 may vary the polish table speed, the wafer carrier down force, the polish time, 
the dimensions or sweep range of the carrier, or the like. Task 96 may also hold one or more 
CMP process parameters fixed (in response to a user input or automatically) to simplify 

30 subsequent processing- 
After task 96 adjusts the initial CMP process parameters, a task 98 causes modeling system 
10 to perform the modeling routine wilh the updated CMP process parameters in place. The 
modeling routine was described above iji connection with FIG. 2. 1 asks 90 and 92 arc repeated 
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to analyze the new simulated CMP result. In this manner, tasks 90, 92, 96, and 98 form a 
processing loop during which the CMP process parameters are optimized in accordance with the 
intended CMP results. Tlie panicular set of CMP process parameters may be saved in database 
34 (see FIG. 1) for future reference or to initialize a subsequent optimization routine. 

.5 Those skilled in the art should appreciate that optimization process 28 (or a modified 

version thereof) may be employed to optimize CMP process parameters related to a theoretical 
CMP system. Such simulations may facilitate the design and development of new CMP systems 
without the cost and labor associated with actual experimentation and prototyping. For example, 
after a particular CMP procedure has been optimized for an existing CMP system « process 28 

1 0 may be performed to var^^ at least one of the optimized CMP process parameters to thereby 
define an updated CMP process parameter set. Then, a second simulated CMP result may be 
obtained using the updated CMP process parameters. In this manner, an operator can efficiently 
simulate and optimize performance of a new CMP system prior to building a prototype. 

It should be noted tlmt optimization process 28 is not limited to the particular optimization 

1 5 protocol described above. Indeed, process 28 may utilize any suitable optimization technique 
to accomplish the same results. For example, process 28 may utilize minimization or 
maximization methodologies to optimize one or more values associated with quantifiable CMP 
results. In the preferred embodiment, process 28 may endeavor to maximize the quality of the 
CMP result by adjusting the CMP process parameters in accordance with the following 

20 technique. A quality measurement (0 is defined as follows: Q - (a)(uniformity) + (]- 
a)(planarity), where the user selects a value for a (the relative importance of uniformity 
compared to planarity) between zero and one; Tlius, process 28 may iieratively adjust the CMP 
process parameters to maximize subject to the particular value of a. 

In the preferred embodiment, modeling system 10 is capable of validating the modeling 

25 parameters associated with a simulated CMP result such that the error between the simulated 
CMP result and an empirical CMP result (using the same CMP process parameters) is 
substantially minimized. Such model validation is desirable to determine whether a given model 
can reproduce experimental data through manipulation of its modeling parameters, to establish 
rules for computing the modeling parameters for new configvirations, and to determine the 

30 reliability of such rules and the accuracy of the simulations. 

FIG. 8 is a flow diagram of model validation process 26 that may be performed by 
modeling system 10. Generally, process 26 enables modeling system 10 to validate its 
simulation results by comparing the simulated results to corresponding empirical results and 
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adjusting the modeling parameters to increase the accuracy of the simulation. Thus, modeling 
system 1 0 may be configured to leverage historical empirical data to improve the performance 
of the CMP simulations. Process 26 may be perfonned for any number of modeling parameters 
and for different CMP process parameters to enable accurate simulations for a wide variety of 

3 different processing environments. 

Model validation process 26 may begin with a task 1 12, which causes modeling system 
10 to retrieve a particular or current CMP simulation. Task 112 may be similar to task 86 
described above in connection with FIG, 7. The CMP simulation retrieved during task 1 12 may 
have been generated by modeling system 10 in accordance with the present invention or, 

10 alternatively, by any suitable modeling methodology. A CMP procedure is conducted during a 
task 114; lliis CMP procedure is perfonned by CMP system 12 in accordance with the CMP 
process parameters designated during the corresponding CMP simulation. Task 1 14 may be 
performed in response to control signals from modeling system 10 or in response to operator 
inputs at CMP system 12. In addition, task 1 14 may be performed at any time before or after 

15 task 112. In an alternate embodiment, task 1 14 may be performed in an interactive manner as 
modeling system produces a substantially real-time CMP simulation. 

Following task 1 14, a task 1 1 6 may be performed to measure an empirical CMP result 
associated with the wafer processed during the CMP procedure. The empirical CMP result is 
preferably stored in database 42 (see FIG, 1) for subsequent use by modeling system 10. In the 

20 preferred embodiment, task 1 16 may be performed by a measurement system incorporated into 
CMP system 12 and/or modeling system 10, or by one or more separate measurement systems. 
CMP system 12 may include any number of in-situ wafer measurement devices to facilitate 
substantially real-time optimization of the CMP process parameters in response to current 
simulation results. Such measurement devices (and other suitable wafer measurement systems) 

25 may be employed for purposes of task 1 16; such systems and devices* which may be known to 
those skilled in the art, are not described in detail herein. 

Task 1 1 6 preferably measures at least the global uniformity of the processed wafer (derived 
from a film thickness profile) and the local planarity of the processed wafer (derived from a local 
feature pattern profile). Thus, task 1 16 obtains an empirical CMP result that contains a wafer- 

30 scale empirical CMP result and a feature-scale empirical CMP result. Thus, the global 
uniformity may be derived by measuring the film or wafer thickness at a number of points on the 
surface of the wafer. In an exemplary embodiment, the thickness profile is measured with a thin- 
film iliickncss measurement system, such as an Optiprobe system. The local planarity may be 
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derived by scanning the surface of the wafer and analyzing the small scale features. For example, 
the preferred embodiment utilizes a surface profile measurement system, which may be 
commercially available as a stand-alone system. Other suitable techniques may be utilized 
during task 1 16 to measure uniformity, planarity, or other characteristics of the wafer, e.g,^ 
5 systems tliat employ reflective or refractive optics or systems that employ micrometer or 
observational techniques. For the sake of compatibility, the measurement points preferably 
correspond to the sampling points used by modeling system 10 to produce the current simulated 
CMP result. 

After task 1 1 6 obtains the requisite amount of empirical CMP data, a task 1 1 8 is performed 

10 to compare the simulated CMP result to the empirical CMP result. In practice, modeling system 
10 may compare the simulated and empirical film thickjiess profiles and the simulated and 
empirical local feature profiles. During task 1 1 8, the wafer-scale empirical CMP result and the 
feature-scale empirical CMP result are compared to their respective simulated counterparts. 
Modeling system 10 may use any suitable technique during task 1 1 8 to compare the empirical 

15 and simulated results. For example, the various sampling points may be analyzed on an 
individual basis or a plurality of sampling points associated with a given measurement may be 
processed in a collective manner. Alternatively, task 118 may employ any number of 
conventional curve fitting techniques, least-squares techniques, or the like. 

In accordance with an exemplary embodiment of modeling system 10, a query task 120 is 

20 performed after the empirical and simulated CMP results have been obtained. Query task 120 
determines whether a simulation error (which may be obtained during task 11 8) is substantially 
minimized. Such a simulation error may be determined for individual sampling points or for a 
collected or averaged quantity associated with a number of sampling points. For example, in the 
preferred embodiment, the difference is mathematically determined by summing the individual 

25 differences between die simulated and empirical results at various data points. Tlie sum of these 
individual differences is dependent upon the values of the modeling parameters used by 
mudeiing system 10. 

If query task 120 determines that the current simulation enor has been minimized, then a 
task 122 is performed lo store the current modeling parameters for use with a subsequent 
30 modeling routine. In other words, the optimized modeling parameter set may be used to conduct 
fliture simulations in a confident manner. The optimized modeling parameteris may be stored in 
database 36 of memory 22 (see FIG. I). After completion of task 122, model validation process 
26 ends. It should be appreciated that process 26 may be performed in conjunction with 

18 



Wo 99/25520 



PCTAJS98/24231 



optimization process 28 in an iterative or combined manner to obtain optimized modeling 
parameters and CMP proce^ss parameters for a given CMP procedure. 

If query task 1 20 determines that the current simulation error is unacceptable, then model 
validation process 26 leads to a task 124. Tosk 124 causes modeling system 10 to adjust at least 
one modeling parameter to obtain an updated modeling parameter set. Modeling system 10 may 
be suitably configured to systematically adjust the modeling parameters in a number of ways and 
in any order. Task 124 may also cause modeling system 10 to hold at least one of the modeling 
parameters fixed during the adjustment procedure to facilitate efficient and speedy optimization. 
The specific modeling parameters to be fixed may be designated by the user prior to model 
validation process 26, during process 26, or by default. It should be noted that, because the 
simulation error is dependent upon the various modeling parameters, process 26 is reduced to a 
multi-variant optimization problem. 

Following task 124, a task 126 preferably causes modeling system 10 to perform a 
subsequent modeling routine to obtain an updated simulated CMP result associated with the 
updated modeling parameters. Task 126 is similar to task 98 described above in connection with 
FIG* 7. After task 126 obtains the new simulated CMP result, model validation process 26 is 
reentered at task 1 18 to re-compare the empirical CMP result with the updated simulated CMP . 
result. Thus, tasks 118, 120, 124, and 126 arc preferably repeated until the error between the 
simulated CMP result and the empirical CMP result has been substantially minimized. In other 
words, the simulated CMP result is altered until a "best fit" relative to the empirical CMP result 
has been obtained. This processing loop causes modeling system 10 to self-optimize its 
modeling parameters such that the simulated CMP result substantially matches the empirical 
CMP result. Although the present invention employs non-linear regression techniques to 
optimize the modeling parameters, any number of suitable metiiodologies may be utilized to 
determine the particular modeling parameters used for a given simulation. 

FIG. 9 is a flow diagram of interpolation formula generation process 30 performed by an 
exemplary embodiment of modeling system 10. Process 30 preferably begins with a task 132, 
during which modeling system 1 0 obtains a plurality of optimized modeling parameter sets, each 
being associated with a specific set oi" CMP process parameters. In other words, each of tiie 
optimized modeling parameter sets is configured for use with the CMP modeling routine to 
obtain a simulated CMP result for a wafer processed during the specified CMP procedure. The 
optimized modeling parameters are preferably obtained during model validation process 26 (see 
FIG. 8) or by an equivalent parameter estimation technique. Task 132 may obtain the modeling 
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parameters from database 36 (see FIG. 1). 

Task 132 leads to a task 134, which causes modeling system 10 to develop a plurality of 
interpolation formulas associated with the optimized modeling parameter sets. The interpolation 
formulas may be utilized by modeling system 10 to predict modeling parameter values for new 
and untested sets of CMP process parameters. Those ^kiUsdjtrt he art w ill appreciate that any 
number of m ethods may be used to form the interpolation rules for the modeling param eters. 
The most straightforward maimer is to simply express each modeling parameter as a linear or 
quadratic function of each of the CMP process parameters. Such a metho dology may be 
acceptable when a large amount of empirical data is accumulated to encompass an adequate 
number of CMP procedures. However, this tactic may not be adequate when a limited amount 
of empirical data exists. 

The preferred embodiment endeavors to generate interpolation formulas associated with 
those CMP variables that dictate performance and/or provide for reliable interpolation rules for 
the modeling parameters. One technique that accomplishes this goal limits the number of 
variable CMP process paramereST:oThe foUowingi 

r= table speed 

C = carrier speed 

F = carrier down force 

Modeling system 10 uses a combination of these parameters for interpolating tlie modeling 

parameters; this combination is dermed as IT^ where: 

^ 

If X* is an objective function of a minimization problem and X" ^2* ^3* ^v) then an 

interpolation rule may be obtained by calculating the different possibie combinations and 
analyzing x~. As an approximation to this approach* the exponents can be restricted to the values 
- 1 , 0, or 1 (with ffj 1^ I). The best overall "fit" is then selected for purposes of the interpolation 

rormOla. X^^^^. 

After the_interpplation formulas are adequately developed in task 134, a task 1 36 may be 
performed such that modeling system 10 receives CMP data associated with a particular CMP 
procedure. As described above, this CMP data may be related to the CMP process parameters 
and/or the intended CMP results {e,g.y uniformity versus planarity balance, die pattern, or the 
like). Task 136 may be pertbmied in respoiisc to user-defined inputs or in accordance with other 
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instructions received by modeling system 10. 

A task J 3 8 is preferably performed in response to the CMP data received in task 136. Task 
138 causes modeling system 1 0 to generate an interpolated modeling parameter set based upon 
the current CMP data. Task 138 utilizes one or more interpolation formulas to generate the 
5 interpolated modeling parameters for the given CMP data. Following task 138, process 30 may 
end. The interpolated modeling parameters may be stored or immediately used to obtain a 
simulated CMP result for the specified CMP procedure. 

As depicted by the ellipses and connecting arrow in FIG, 9, a task 140 (which may be 
pertbmied at a later time or in connection with a separate process) may be perfonned to produce 

1 0 a plurality of CMP process parameters in response to an intended CMP result, where the CMP 
process parameters are obtained through the optimization procedures described above, which 
involve the consulting of tJie interpolation formulas. Thus, process 30 (or a related process) may 
be performed by modeling system 10 to obtain suggested CMP process parameters for an 
untested CMP procedure if an intended CMP result is known beforehand. 

15 The CMP process parameters obtained during task 140 may be subsequently (or 

substantially concurrently) applied to CMP system 12, as depicted in a task 142 related to 
process 30. Task 142 may cause modeling system 10 to communicate the CMP process 
parameters to CMP controller 32 within CMP system 12 (sec FIG. 1). Alternatively, an operator 
may make a record of the CMP process parameters and adjust CMP system 12 in an appropriate 

20 mariner to effect the desired settings. Following task 142, a task 144 causes CMP system 12 to 
perform a CMP procedure in accordance with the CMP process parameters produced during task 
1 40. Task 144 is similar to task 1 14 described above in connection with FIG. 8. Following task 
144, process 30 ends. As described above, the empirical CMP results obtained during task 144 
may be used during model validation process 26. 

25 In summary, the present invention provides an improved CMP modeling system that is 

capable of providing both wafer-scale simulations and feature-scale simulations. The CMP 
modeling system can generate an optimized set of CMP process parameters based on a specified 
balance between wafer uniformity and planarity. In addition, the modeling system computes its 
modeling coefficients in accordance with empirically determined CMP results to thereby reduce 

30 simulation errors. The modeling system employs interpolation techniques to hs modeling 
parameters to effectively simulate CMP results for which little or no empirical data exists, and 
it can process simulation data collected for an existing CMP system to assist in the design of a 
new CMP system having different physical and processing characteristics than the existing CMP 
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system. 

The present invention has been described above with reference to preferred exemplary 
embodiments. However, those skilled in the art will recognize that changes and modifications 
may be made to the preferred embodiment without departing from the scope of the present 
5 invention. For example, the present invention is not limited to the particular modeling, 
optimization, or interpolation techniques described herein. In addition, the various hardware 
components may differ than that shown and described herein and the various processes need not 
be performed in the precise manner described herein. These and other changes or modifications 
are intended to be included within the scope of the present invention, as expressed in the 
10 following claims. 
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